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We study the corrections to e'/e in the minimal supersymmetric model at the leading order 
in QCD and QED. Supersymmetry can increase the standard model prediction for e'/e by 
at most 40% for m t = 174 GeV, an enhancement which is indistinguishable from the present 
theoretical uncertainties. The most conspicuous effect of supersymmetry is a strong depletion 
of e'/e. For certain choices of supersymmetric parameters, vanishing and even small negative 
values of e'/e can be obtained for the top quark in the CDF range. 
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1 Introduction 

CP violation in the kaon system has always been a fertile field in which to test theories beyond 
the Standard Model (SM). At present the situation for e'/e is particularly interesting. The 
experimental results from NA31 at Cern |lj 

Re e'/e = 23 ± 6.5 x 1(T 4 (1) 

and from E731 at Fermilab || 

Re e'/e = 7.4 ± 6.0 x 10" 4 (2) 

are incompatible at the l-a level. On the theoretical side, there has been great progress 
in reducing the uncertainties in the prediction for e'/e. The full next-to-leading QCD per- 
turbative corrections have been computed by two groups |3], |J] and developments in lattice 
calculations should soon provide reliable estimates of the hadronic matrix elements. 

Given the inconclusive experimental situation and the notable refinement in the SM 
calculation for e'/e, the time is right for detailed investigations of e'/e in theories beyond 
the SM. If the top quark is indeed as heavy as implied by the CDF measurements || , the 
next-to-leading SM prediction for e'/e sits comfortably in the E731 range, disfavoring the 
larger NA31 result. While we await improved experimental statistics, which will resolve 
the conflict between eq. ([!]) and eq. (|2|), it is now interesting to see if extensions of the SM 
can substantially modify the prediction for e'/e and possibly account for the large value 
suggested by the NA31 measurement. Furthermore, since the SM prediction for e'/e suffers 
a strong accidental cancellation among the different contributions for the top quark in the 
range of interest, it is conceivable that even modest new- physics effects will easily stand out 
and sizeably modify the final result. 

A notable example of study of e'/e in theories beyond the SM is contained in ref . || . There, 
in the context of the two-Higgs doublet model, a systematic analysis of e'/e including leading 
QCD logarithms has been carried out. It was found that the charged Higgs has the effect 
of reducing e'/e with respect to the SM prediction. This reduction occurs partly because 
of a positive charged-Higgs contribution to e, which suppresses the Cabibbo-Kobayashi- 
Maskawa (CKM) CP violating parameter sin S, and partly because of a new contribution to 
electroweak penguin and box diagrams, which enhances the degree of cancellation with the 
strong penguin already present in the SM for a heavy top quark. 

The aim of this paper is to generalize the analysis of ref. || to the case of the minimal 
supersymmetric model. Although several analyses of e'/e in the context of supersymmetry 



are already present in the literature J?) , a complete study including leading QCD logarithms 
has not been performed. In order to attempt such a study, it is necessary to make precise 
assumptions regarding the structure of CP violation and flavor- changing neutral currents 
(FCNC) in the supersymmetric model. Here we choose to work in a minimal version of the 
model, described in more detail in sect. 2, where both CP violation and FCNC are completely 
determined by the usual CKM matrix elements. We have several reasons to do this: 

i) The corrections to e'/e considered here are fairly generic to all supersymmetric models. 
Non-minimal models may contain new CP-violating phases and/or tree-level FCNC, thereby 
generating extra contributions to e'/e. These contributions are however strongly model- 
dependent. 

ii) The minimal version of the supersymmetric model considered here is very predictive, 
since CP violation and FCNC are completely determined in terms of only the known CKM 
angles and phases. 

in) The effective Hamiltonian below the weak scale can be written in terms of the same 
set of operators used in the case of the SM. This greatly simplifies the analysis, since the 
anomalous dimension matrix is then completely known at the leading || and next-to-leading 
U order. If, for instance, flavor-changing gluino-mediated interactions were included, new 
operators with different chiral structure would appear and a calculation of new elements of 
the anomalous-dimension matrix would be required. 

Our paper is organized as follows. In sect. 2 we describe the version of the supersymmetric 
model under investigation and we establish our notation. In sects. 3 and 4 we give the 
formulae for the supersymmetric corrections to e and e'/e. The results of our numerical 
investigation are presented in sect. 5. Finally, in sect. 6 we summarize our results. 

2 Supersymmetric model with minimal CP and flavor 
violation 

In general, supersymmetric models have potential sources of new CP violation because of 
the presence of unremovable phases in the supersymmetry-breaking terms. However, these 
phases are strongly constrained by measurements on the electric dipole moment of the neu- 
tron and must be small, typically less than 10~ 2 — 10~ 3 Q. Although there is no compelling 
theoretical argument which suggests that they are exactly zero, we will make the simplifying 
(and often used) assumption that the only CP violation resides in the CKM matrix. 

It is also well known that FCNC in supersymmetric models can arise at the tree level, 



since the quark- and squark-mass matrices are diagonalized by different field transformations 
fllOfl . Once the heavy supersymmetric particles are integrated out, one is left with FCNC 
involving ordinary quarks, which are induced by one-loop diagrams with squark and gluino 
exchange. Their effects are potentially large, since they are 0(a?), as opposed to the stan- 
dard model FCNC contributions O(a^). The flavor-changing quark-squark-gluino couplings 
depend however on the details of the supersymmetry-breaking sector, which is still the least- 
understood part of the theory It is customary to make the simplifying assumption that all 
supersymmetry-breaking terms are flavor independent at some grand-unification scale, close 
to the Planck mass. This is the case if, for instance, the Kahler metric of the underling 
supergravity theory is flat. In this case flavor- changing quark-squark-gluino couplings are 
induced by renormalization effects, but their influence in kaon physics is negligible because 



of the strong limits on gluino and squark masses from hadron colliders [11 



There is growing criticism of this point of view ||12| , since supergravity theories derived 
from superstrings do not seem to have flat Kahler metrics and do not seem to have flavor- 
independent supersymmetry-breaking terms. However, if the supersymmetry-breaking terms 
have a completely general structure in flavor space, they lead to corrections to K° — K° 
mixing larger than that allowed by experimental constraints. In this case there is need for a 
suppression of the flavor- changing quark-squark-gluino coupling either by postulating some 
form of flavor symmetry at the unification scale or by invoking some dynamical mechanism 
to make squarks more degenerate in mass (e.g. the running from unification to weak scale 
in a gaugino-dominated supersymmetry-breaking scenario). 

In our analysis we will not be concerned about the fundamental mechanism which sup- 
presses FCNC, but we will simply postulate, as we have done for CP violation, that, at the 
weak scale, FCNC are absent at tree level. As discussed in the introduction, this hypothesis 
of minimal CP and flavor violation enhances the predictivity of the model and largely sim- 
plifies the analysis of QCD effects, since the operator basis of the effective Hamiltonian is the 
same as in the SM. Let us now briefly present the theoretical framework in which we work, 
the supersymmetric standard model with minimal CP and flavor violation, and establish our 
notation. 

By making a transformation on superfields, we choose a basis in which the up-quark mass 
matrix m u is real and diagonal and the down-quark mass matrix is Vrrid, where rrid is real 
and diagonal and V is the CKM matrix. In this basis, the 6x6 up- and down-squark mass 
matrices are: 

Ml = ( a Y ml -M\). (3) 

V A um u m z UR + m z u J 



T-,2 



m 



VmlV^ A r ,Vm n 



M 2 d = V T T'""\]- (4) 

d \ A d m d V^ rh dR +m d ) K) 

Our minimal CP- and flavor- violation hypothesis states that, at the weak scale, the supersym- 
metry-breaking terms rh UL R , vn dh R and A U)d are CP-conserving and flavor- independent, i.e. 
real and proportional to the 3x3 identity matrix. If this hypothesis holds, there are no 
tree-level flavor-violating quark-squark-gluino couplings at the weak scale. Because of the 
strong constraints imposed by the real part of K° — K° mixing, we know that this hypothesis 
should be approximately satisfied, unless squarks or gluinos are very heavy. 
The mass matrices eqs. (§)-© can be diagonalized according to 

rn iC^rri j- ( m u Li Jrm u m . 2 i A u m Ui m +m 2 A u m u \ , . 



'i " "■"■" ""i 



™2 i ™2 a ™ ™2 i ~ 2 



lll \i+ m d Ri j_ _ 2 ^ Atm* ™.d Li + ™>d Ri ^ _ 2 __ A d m dt 
with % = 1, 2, 3 and where 

-s„ a,)' '" V-&V* o. 



TMlTl = d.ag -*^* + mjj, + iS, ^Y^ + < " j£S I , (6) 



T ^ = o T ,+ ^ , ( 7 ) 



'u 



C u4 = diag (cos 6»"' d ) , S^ = diag (sin 6^' d ) , (8) 



9 4 rn U ' d 

tan2C' d = , 2 M J 2 . (9) 

" b u Li ,d Ll " L u m ,d Ri 

The chargino mass matrix is diagonalized by two orthogonal 2x2 matrices U and V, 
according to: 

u( M r m ^ si ^)v -=(*» M, (10) 

\m w v2cos/? [i j \ m X2 / 

where tan /3 is the ratio of the two Higgs vacuum expectation values, M is the weak gaugino 
mass, and /i is the Higgs superfield mixing parameter, taken here to be real. 

In order to compute the relevant one-loop diagrams, we need the interaction of charginos 
(Xjj J = 1) 2) with down quarks (d) and up squarks (ilk, k = 1, 2) in terms of mass eigen- 
states: 

C x = gdV^ (z* 1 -^ + yi*I±2*) ^u k + h.c, (11) 




Z 3k = _^£ | | (12) 

V2m w cos (3 ° ' ■- "" ' 

Y ' k = 7&P V » l "<-> "^ """-'■ (13) 

where flavor indices are understood. The CP- and flavor-violation properties of the interac- 
tion Lagrangian in eq. (|Tl|) are determined by the familiar CKM matrix. This property is 
satisfied also by the charged-Higgs interactions with quarks: 

C H = -jS—H+u (J—m u V^P^ + tan (3VmJ-^- J d + h.c. (14) 

\/2m w \tanfj 2 2 ) 

We conclude this section by summarizing the free parameters necessary to specify the 
supersymmetric model, in addition to those of the SM. Whenever convenient, we choose to 
redefine the original parameters and work with physical particle masses as input. For the 
first two generations of up-squarks m u is negligible with respect to fh UL R and we will take 
these squarks to be degenerate in mass, i.e. m UL ~ fh UR = m. This is merely a simplifying 
assumption which will not affect our conclusions. The mixing between the two stops cannot 
be neglected. We choose as input parameters the lightest stop mass, mj, and the stop mixing 
angle 9 (which can be chosen in the range between and tt, without loss of generality). The 
mass of the heavier stop (rrif) is then given by: 

rrVf = 2m + 2m^ — m~. (15) 

The chargino mass matrix is defined by three parameters: we choose them to be the lightest 
chargino mass (m x ), the ratio of Higgs vacuum expectation values (tan/?), and the weak 
gaugino mass (M). The last free parameter is the charged-Higgs mass (Mjj+), while the 
relevant H ± couplings are completely determined once the value of tan (3 is fixed. Therefore 
the model requires seven input parameters: m, mj, 8, m x , tan/3, M, Mh+- Notice that 
some of these seven free parameters may become related with one other if additional assump- 
tions are made (specific relations at the GUT scale, constraints from electroweak symmetry 
breaking, etc.). We prefer to treat them as independent variables in order to describe a 
larger class of models. 



3 The parameter e 



The parameter e gives a measure of the indirect CP-violation in kaon decays. It is defined 

as 

e wr/4 



• - -, "ImM 12 , (16) 

y/IArriK 

where M 12 is the imaginary part of the off-diagonal element of the neutral kaon system 

(K°, K°) mass matrix, and Am^ is the Kl — K$ mass difference. 

By evaluating the relevant AS = 2 box diagrams in figure la-b, we obtain 



in/4 G F 2 m I 

12V2n 2 K A mi 



-B K m 2 w Y,^ 2 c)v?S ii 



[c,c) + 



(lm\ 2 t )4 l) S^\t,t) + 2(ImA c A t )^ ) S®{c,t) 



(17) 



where A» = V* d Vi S , f K = 161 MeV is the kaon decay constant, and Bk is the non-perturbative 
parameter which defines the normalization of the relevant hadronic matrix element in units 
of the vacuum-insertion value. The coefficients rji include the QCD corrections to the box 
diagrams and are given in table 1. The values of r\i for the box diagrams involving squarks 
and charginos are all equal to the SM value of 772 because we are assuming that the super- 
symmetric particles are integrated out at the scale of the W boson. 

The functions S^ 1 ' (c, t) result from calculation of the box diagrams with exchange of 
quarks and two W bosons (jyW 7 ), two charged Higgs bosons (HH), one W and one charged 
Higgs (HW), and with exchange of squarks and charginos (x)'- 



qWW 

S HH i 



c,t) 
c,t) 



S{x 



LoiX 



4tan 4 /3 



CWl ■Etwi 



s 



HW (c,t) 

S x (c,t) 
f(c,t) 



tan 2 /3 



- 1j2{X cw , Xtwy %I 



L'X \&CW 1 %tWl %HW ) 



f(u,u)-f(u,c)-f(u,t) + f(c,t) 



E 

h,k=X,2 



X 



W X j 



^ichXit^jCfXjik L 3 \%i k 



Xji Xc hXji X XiXi)- 



(18) 



We have denoted x a b = m^/m 2 . for generic indices a, b, and have collected the functions S 
and £1,2,3 i n the appendix. The expressions for the charged-Higgs contributions to short- 
distance effects given in this and the next section agree with ref. 0, and the expressions for 
the chargino contributions agree with the existent literature whenever results are available 
(see, in particular, ref.||13|). 
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For a given set of supersymmetric parameters, one can compare the theoretical prediction 
for e in eq. fllTD with the experimental result and then determine S, the CP- violating phase in 
the "standard" representation [[I i] of the CKM matrix. In general one obtains two solutions 
for S but, in most of the parameter space, one of them is ruled out by the experimental data 
on B° — B° mixing and on the b — > sj branching ratio. 

The B° — B° mixing is determined by the box diagrams in figure lc-d in a manner 
analogous to e: 

x d = ^ = r B ^ VQCD m B B B f B m 2 w \V td \ 2 j:S®(t,t). (19) 

We take m# = 5.28 GeV, the perturbative QCD correction parameter t]qcd = 0.84, and the 
experimental determination from ARGUS |I5| and CLEO [|T(J for x d = 0.70 ± 0.13. 

The expression for the decay rate of b — > 37 in supersymmetry can be found in ref. 
[U| (or, in a notation very similar to the one employed here, in ref. |17|| ) and will not be 



repeated here. In our analysis we will impose the constraint on the inclusive branching ratio 
BR(b — > sj) < 5.4 x 10~ 4 obtained by CLEO fI8f . The numerical results will be presented 
in sect. 5. 



4 The parameter e' 



The parameter e' gives a measure of the direct CP- violation in kaon decays. It is defined as 
follows: 

(20) 



e = 


--> 


-Q) 


e<*, 


where 








i = 


lmA 


UJ = 


ReA 2 



^ 1 ImA 2 Ti „ „ a , s 

p , P , fi = -r^' <p=~+$2-So^j. 21 

ReA ReA uj IthAq 2 4 

Aj and 5j are the amplitudes and final state interaction phases for the decay K — > (7171)1, 
where I denotes the isospin of the final pion state. Experimentally one has 

ReA Q = 3.3 x l(T 7 GeV, uj ~ 1/22. (22) 

The procedure to relate the amplitudes in eq. (|2"DD to the QCD-improved AS" = 1 effec- 
tive Hamiltonian is by now standard flT9| , pl|] We will follow the method illustrated in ref. 
pl|, briefly review its results, and show where modifications due to the supersymmetric 



contributions are necessary. 



In the SM, after integrating out the heavy degrees of freedom at the scale of the W boson, 
the tree-level AS = 1 Hamiltonian becomes: 



Gr p 



where 



V2 



v ud v* 



Hffr 1 = -^V ud V: s (1 - r) {Ql - Q c 2 ) + r Q« - Q 



Ql = {sq) V -A (? d ) 



V-A VH^lV-A ' 



q — u, c,t. 



(23) 



(24) 



With QCD corrections taken into account, the effective Hamiltonian at an energy scale 
fj, below the charm mass is obtained by replacing in eq. ( p3|) the operators Ql with the 
renormalized four-fermion operators Qi{[J>) as follows: 



10 



10 



(Q u 2 - Ql) -+ £*(aOO«(a«), (Q u 2 - Ql) - £«MG<M- (25) 

»=1 i=l 

In eq. fl25|) , fj(/f) and £j(/i) correspond to the relevant Wilson coefficients and the complete 
basis used for the operators Qi is 



Qi = (sd) v ^ A (uu) v _ A , 

Qi = {sd) V -A J2 (QQ)v-a, 

q=u,d,s 

Q5 = (sd)v„ A £ (qq)v+A, 

q=u,d,s 

3 

~ q=u,d,s 



Q2 = {su) v _ A (ud) v _ A , 
Qi = £ (sq)v-A(qd)v-A, 

q=u,d,s 

-8 £ (s L q R )(q R d L ), 

q=u,d,s 

Q 8 = - 12 £ e q (s L q R )(q R d L ), 



Qt 



q=u,d,s 



Q9 = -(sd) v _ A £ e q (qq) v _ A , 

~" q=u,d,s 



Qw = j £ e g (sg)i/_ A (grf) y _ A , (26) 

where e g is the electric charge of the quark g, (V ± A) refer to 7 M (1 ± 75), and qi jR = 
1/2(1 ±75)5. We neglect here magnetic-moment operators [22], since their contribution is 
unimportant for our purposes. 



In general we can therefore write the effective AS" = 1 Hamiltonian at the scale \x as 
G 10 



H. 



where 



v2 i=1 



(27) 



CiQi) = Zi(ji) + Tyi(p)] yiQi) = Vi(p) - z t (fi). 



(21 



The Wilson coefficients C{ are chosen to match the short- distance contribution at the ap 
propriate scale where heavy particles are integrated out, and then evolved using the renor 
malization group equation: 

d „, x d 



^ + fo)f g -?&><*> 



i 2 
m w „2 



Cl-f,^,al=0, (29) 



where C stands for a ten-dimensional vector, f3(g) is the QCD beta function, and a is the 
electromagnetic coupling constant (the running of a is neglected). In this paper we use the 
anomalous-dimension matrix 7 T (g 2 , a) at the leading order in QCD and QED as given in ref. 



n[|, although it is now known at the next-to- leading order ||. This is perfectly adequate, in 
view of the large intrinsic uncertainty due to the ignorance of the supersymmetric parameters 
and to the inescapable model-dependence. Our aim here is just to illustrate the trend of the 
supersymmetric corrections to e'/e. 

The modifications caused by supersymmetry appear only in the boundary conditions 
of the Wilson coefficients, which are computed through the appropriate one-loop Feynman 
diagrams in figure 2. In our analysis we will impose the boundary conditions at the scale 
\i = m w , although strictly speaking they should apply to the energy scale at which the heavy 
particles are integrated out. The threshold corrections, originating from this mismatch of 
energy scales and certainly present in any realistic theory with a non-degenerate mass- 
spectrum, are numerically not very significant, since the running of a s in the region above 
m w is not very steep. In the HV renormalization scheme, the boundary conditions for the 
Wilson coefficients in the minimal supersymmetric model are given by: 

/ x a s (m w ) 
vx{m w ) = — (14 -Bg) 

107T 



v 3 (m w ) - 



a l f B W | B(U) | A < m ^ E 



67rsin 2 # w \ 2 I 24tt 



/ v a s (m w ) 

Vi{m w ) = — E 

f v a s {m w ) 

V(i{m w ) = — E 

v 7 (m w ) = ^(AC + D) 

07T 

v&{m w ) = 



V9(m w , 



a 

6tt 



AC + D + — 2 (- B{d) + B{u) ~ 4 C 

sin 9 W v 



vio(m w ) = 



(30) 



and 



zi(m w ) = vi(m w ), z 2 (m w ) = v 2 (m w ), Zi(m w ) = 0, i — 3, 10. 



(31) 



The functions which include the contributions from photon-penguins (D), Z-penguins (C), 
gluon-penguins (E), boxes with external down quarks (B^) and up quarks (B^), and 
gluino-mediated boxes (B^) are given by: 



D = D SM (x tw ) + T7 -^D H (x tH ) + Y, [ Y j 2 i k x wi k D x( x Xj i k ) ~ ' f _> r 



tan 2 (3 



fe=l,2 



(32) 



x^ 



tan^p 



«,j=l,2 

h,k=l,2 



C H (x tH )+ Y, YAY tih \-5 tJ A hk C£\x ihX] ,x ik 



'*Xj> 



+ 5, 



/iA: 



UixUjxCMix^x^) + VnVj! — logm| - Cj'^.^^- 



£ ^c 



£ = E SM (x tw ) + —i r -rE H (xt H ) + Yl [ Y k x wi k E x( x Xj i k ) - (* ~* 5 ) 



tan 2 /3 



3=1,2 
fc=l,2 



(33) 
(34) 



B (d) = -2#5m(^) + - E ******* 



»,j=l,2 
h,fc=l,2 



*j*fc *#fc -°X v X *feXj ' X «hXj ' X XiXj - 



t-^C 



(35) 



5 (M) = 8B SM (x tw ) + - J2 UjiUii^wxiVx^lyjiJii^Kx^^Xi ,x XiXj ) 



i,3=l,2 



(f^a) 



(36) 



-D — 2^ ^/ifc^ivg \^ il "^ »lj X y X Xigi X Ukgi X u h g) 

i=l,2 

h,k=l,2 

+ 2^^V,U,B^(t t- t- \ 



m n 



(37) 
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where Ahk = 1 for h = k and Ahk = —1 for h ^ k. Finally, the explicit form of the functions 
resulting from loop integration is given in the appendix. 

Notice that the supersymmetric box diagrams contribute to the operators Q3 and Q 9 via 
two different functions B^ and B^ d \ as opposed to the SM case where a single function 
(Bsm) appears. The function B^> accounts for the contribution of the box diagrams with 
a gluino, a chargino, and two squarks as internal lines. These diagrams exist also in the 
limit of flavor-conserving gluino vertices, in which we are working. In order not to further 
increase the number of free parameters, we have chosen to relate the gluino mass (m§) to 
the M parameter defined in sect. 2 through the GUT relation: 

a s (m w )sm9 w M 

m~ g = . (38) 

a{m w ) 

This assumption is however inessential, because the gluino-box contribution to e'/e is com- 
pletely negligible. 

Using the effective Hamiltonian H^ s l =1 of eq. (|27|) , we recast eq. (|20|) in the form 

- = ImA *9i^r^0») (QM)o (1 - fi ) (39) 

where £l v +rf is the contribution from 7r — rj — t]' mixing, here taken to be £l v +rf — 0.25 ±0.05. 
Following a general convention, we have factored out the dominant contribution from the 
operator Q 6 and call Q the correction due to the other operators. 

In our numerical analysis, we set the renormalization scale ft at 1 GeV. In order to 
evaluate eq. (|39|) , we need the relevant hadronic matrix elements at that scale: 

(Q i ) I = ((irirhlQilK) 7 = 0,2, (41) 

where I is the total isospin of the final pion-pion state. In the future, lattice calculations 
should provide the most reliable evaluation of hadronic matrix elements. For the moment 
we prefer to follow the approach of ref . . 

The strategy is to compute the matrix elements (Qi)j at the m c scale and then, using 
the QCD evolution, to obtain the matrix elements at the scale \x < m c . The advantages of 
this method are explained in detail in ref. |J. Starting at ji = m c one finds ||: 

(Qi) = -\XB?'*>, (Q 2 } = 5 -XB^\ (Q 3 ) = \XB^\ 

11 



(Q*) 

(Qe) 

(Q 7 ) 

(Qs)o 

(Q 9 )o 



(Q 3 ) + (Q 2 ) - (Qi) , (Q*) 



1 „n /9i i-^- 



-B 



(1/2) 



Qi 



m. 



m s (/j,) + m d (ji) 

5<55> <* + i>-£ 



' /t r (1/2) 

o(l/2) 
n 7 5 



5 



(V2) 



5 <Qi> - \ (Qs) , (Qia) = (Q2} + \(Qi)o-\(Q3}o 



(42) 



(Qi) 2 

(Qr) 2 

(Qs) 2 
(Q 9 ) 2 



((h), = A ~fxB[^ 

k '7V\ X 



6v/2 x /o y/2 

\ to) + ^X 

2v /2\^ D /o 6 



(Qio) 2 — 2 (^1)2 ' 



5 



(ft) 5 



(3/2) 

7 ) 



5. 



(3/2) 



where 



k 

X 



A 2 

x 



m 2 K — m 2 



4.55 



Qe 



]J^U{<-0 3.71 x K)- 2 G(A^ 



B, 



(1/2) 



z = 3,..., 6, 



(43) 



(44) 



/-1 /q\ fS/2^ 

and 5,- , -Bj stand respectively for the AI = 1/2, 3/2 transitions. 

These matrix elements are obtained using the vacuum-insertion approximation and soft- 
pion theorems, with all of the non-perturbative information included in the Bi parameters. 
We will use the following values at the scale /1 = m c M: 



B\ 



(1/2) 



5.2, 4 1/2) 



B\ 



(3/2) 



0.55, B 



(3/2) 



5.8 ±1.1, B, 
= 1±0.2. 



(1/2) 



1±0.2, 



(45) 



The remaining Bi parameters play only a minor role in the analysis of e'/e and we will set 
them equal to 1. 
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5 Results 



It is well known [^, 21], [§, |J that the SM value of e'/e is typically of order 1CT 4 for 
m t = 150 — 190 GeV, decreases for increasing m t , and becomes zero for m t = 200 — 220 GeV. 
This accidental cancellation occurs because the parameter Q in eq. (|39|) approaches 1 as the 
top quark mass reaches these large values. 

There are two main sources of uncertainty in the SM prediction for e'/e. One affects 
the Wilson coefficients and comes from uncertainties in SM parameters such as the CKM 
angles, m t , m s , and Aqcd- As measurements on these parameters become more precise, the 
calculations of the Wilson coefficients stand on firmer ground. We summarize in table 2 the 
SM parameter values (and their relative errors) chosen for our numerical analysis. 

The second source of uncertainty stems from the hadronic matrix elements, i.e. from the 
Bi parameters discussed in sect. 4. At present lattice calculations seem to be the only tool 
capable of significantly reducing errors. 

Our goal in this paper is to study the effect of supersymmetry on e'/e. Because of the 
cancellation mentioned above, the SM value of e'/e is typically smaller than its separate 
contributions. In this situation, new-physics effects have the chance to emerge and increase 
significantly the total result. 

In our analysis we will vary in the experimentally-allowed region the seven free parameters 
of the minimal supersymmetric model defined in sect. 2. The parameters are subject to the 
constraints that all charged supersymmetric particle masses be heavier than m^/2 and that 



the common squark mass satisfy the CDF bound m > 126 GeV JTTJ. We also use the tree 



level mass bound Mh+ > m w obeyed by the supersymmetric charged Higgs boson and the 
bound tan/? > 1 implied by radiative electroweak symmetry breaking. 

Since we are interested here in studying the trend of the supersymmetric corrections, we 
start by fixing the SM input parameters to the central values given in table 2. As can be 
seen from eq. (|39|) , there are three sources of new-physics corrections to e'/e: the CKM phase 
(ImAi), the strong penguin (y 6 ), and the rest of the electroweak contribution (Q). Let us 
first consider the CKM phase 5, which is determined by the e parameter, eq. fll7|). 

By scanning over the supersymmetric parameter space, we obtain the allowed region of 5 
shown as a function of m SUS y in fig. 3a (for M H + = 80 GeV) and in fig. 3b (for M H + = oo), for 
m t = 174 GeV. We have defined msusy as the mass of the lightest charged supersymmetric 
particle (either the chargino or the stop). With this definition, it is easy to see from fig. 3 
how improvements on the experimental limits on supersymmetric particle masses (e.g. from 
LEP 200) will affect the allowed range of 5. The two SM solutions for cos 5 correspond to 
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the minimum values allowed by the supersymmetric bands in fig. 3b. In other words, the 
presence of supersymmetry has the effect of reducing the value of sin 5, both for the solution 
in the first quadrant (cos 5 > 0) and in the second quadrant (cos 5 < 0); the lighter the 
charged Higgs, the stronger the reduction. Since e'/e is proportional to sin 5, this has the 
effect of decreasing e'/e, as discussed below. If the constraints from B° — B° mixing and 
b — » sj are imposed, only the solution in the first quadrant survives, as shown in fig. 4. 
This is true only because we have fixed the SM parameters at their central values. If all 
uncertainties are properly taken into account, some solutions for S in the second quadrant 
can still be found. In general however a heavy top disfavors solutions with cos S < 0, unless 
Jb is small. 

The impact of supersymmetry on the prediction for e'/e is illustrated in fig. 5. This shows 
the allowed range of e'/e as a function of m SU sy f° r M H + = 80 GeV (fig. 5a) and M H + = oo 
(fig. 5b) with m t = 174 GeV and with the parameters subject to the constraints from B° — B° 
mixing and b — > sj. In fig. 5a the strong reduction of e'/e in the presence of a relatively light 
charged Higgs is apparent. This is partly caused by the decrease in sin S mentioned above. 
However the main reason for such small values of e'/e is that Q can approach 1, allowing 
a complete cancellation for values of m t smaller than that found in the SM. As shown in 
fig. 5a, e'/e can even reach negative values if M#+ is small enough. The charged-Higgs effect 
on e'/e has already been studied by the authors of ref.|| and we confirm here their results. 

As discussed above, the chargino contribution has always the effect to lower sin 5. How- 
ever, depending on the choice of the supersymmetric parameters, the chargino contribution 
to Q can have either sign. As a consequence, e'/e in supersymmetry can also be enhanced 
with respect to the SM value, as seen in fig. 5. This enhancement is numerically rather small, 
at most 40%. For comparison, the SM value of e'/e, which can be read in fig. 5b in the limit 
m susy - ► oo, corresponds to the leading-order calculation. We recall that next-to-leading 
effects tend to suppress even further the value of e'/e 0, [| . 

In tables 3 and 4 we report the values of the three different contributions sin 5, Q, and 
2/6 (Qejo corresponding to the minimum and maximum values of e'/e. Notice that the main 
effect of supersymmetry resides in sin 5 and Q, while y^ (Qe) is never reduced by more than 
5%. Finally fig. 6 shows how the prediction for e'/e is affected as the constraints from B° — B° 
mixing and b — > s'j are removed. 

It is interesting to know for which values of the supersymmetric parameters e'/e reaches 
its minimum and maximum. The charged-Higgs contribution is maximized for the minimum 
allowed values of M H + and tan/5, respectively equal to m w and 1. Notice that all charged- 
Higgs contributions are proportional to l/tan 2 /5 and therefore rapidly decrease as tan/3 
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increases. The chargino contribution is typically maximized when the chargino and the 
lightest stop masses are small (i.e. of the order of wisusy) an d fh — > oo. In this limit, flavor 
symmetry is indeed maximally broken. Also the chargino contribution minimizes e'/e for 
tan/? = 1 and maximizes it for tan/3 in the range 3-5. 

Next we want to take into account the experimental and theoretical errors of the SM input 
parameters. We have used the following approach. We first fix the SM parameters to their 
central values and, for given m SU sy and M H +, we compute the supersymmetric parameters 
which respectively minimize and maximize the contribution to e'/e. Then, keeping these 
values of the supersymmetric parameters fixed, we generate a large number of configurations 
using gaussian and flat distributions for the SM inputs with respectively experimental and 
theoretical errors. We obtain an event distribution for e'/e, from which we can compute 
the average value and standard deviation. Notice that by following this procedure we are 
averaging over solutions with positive and negative cos 5, both of which are allowed when 
the uncertainties on the SM inputs are taken into account. 

The results are shown in tables 5-7, for three values of the top quark mass, m t = 
150, 174, 190 GeV, imposing the constraints from B° — B° mixing and b — > sj. The effect of 
relaxing these constraints is illustrated in table 8. We observe that the enhancement of e'/e in 
supersymmetry is indistinguishable from the SM result within a standard deviation, even in 
the most favorable case m t = 190 GeV. On the other hand, the depletion of e'/e is substantial 
and statistically significant. Unfortunately the experimental sensitivity is presently not 
sufficient to identify such small values of e'/e. 

6 Conclusions 

We have studied the corrections to e'/e in the supersymmetric model with minimal CP and 
flavor violation at the leading order in QCD and QED. Our results are the following. 

Supersymmetry can enhance the SM prediction for e'/e by at most 40% for m t = 174 
GeV and up to 60% for m t = 190 GeV. This enhancement is of the same order of magnitude 
as the SM prediction within one standard deviation and therefore it is not experimentally 
distinguishable. In this respect, minimal supersymmetry cannot explain a value of e'/e as 
large as the one suggested by the NA31 measurement 0. The enhancement is caused by a 
reduction of the electroweak factor Q and it is attained for chargino and stop masses just 
beyond the LEP limit, with all the other squarks and the charged Higgs considerably heavier. 

The most conspicuous effect of minimal supersymmetry is however a strong depletion of 
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e'/e. For certain choices of supersymmetric parameters, vanishing and even small negative 
values of e'/e can be obtained for the top quark in the CDF range ||, m t = 150 — 190 
GeV. Unfortunately sensitivity to values of e'/e below 10 -4 represents experimentally a very 
challenging proposition. This depletion is caused by a reduction of the value of sin 5 extracted 
from e and by an increase in the electroweak correction Q, which can become equal to 1 or 
larger. From the short- distance point of view, the dominant corrections at the origin of the 
e'/e depletion come from a light charged Higgs, and from charginos and stops just beyond 
the LEP limit. 
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Appendix A 



In this appendix we give the expressions for the functions which enter the penguin and box 

diagrams. 

Box(A S=2 ) 



S(x,y) 

Lx(x,y,z) 

L 2 (x,y,z) 

L 3 (x,y,z) 
F(x,y,z) 



xy 



1 3 1 
4 + 21-x 



4 (1 — x) 2 \ x — y 
xy [F(x, y, z) + F(y, z, x) + F(z, x, y)} 

xy [xF(x, y, z) + yF(y, z, x) + zF(z, x, y)] 

— L 2 (x,y,z) 

xy 

xlogx 



logx 3 

+ {x <-> y) - - 



4(l-x)(l-y) 



(x — l)(x — y)(x — z) 



Box(A S=l) 



B S m(x) 



B d (x,y,z) 



x 



+ 



X 



4(x-l) 4(x-l) 2 

L x (x,y } z) 

xy 

L 3 (x } y } z) 



log 2 



16 



7-pengum 


D S m{x) 


D h {x) 


D x (x) 



x 2 25 - 19x) -3x 4 + 30a; 3 - 54a; 2 + 32a; - 8 , 

— r^ + 1 -log a; 

36 (a; -l) 3 18 (a; -l) 4 & 

x (47a; 2 - 79a; + 38) a; (-3a; 2 + 6a; - 4) 



108 (a; -If 18 (a; - l) 4 

-43a; 2 + 101a; - 52) (2a; 3 - 9a; + 6) 



log a; 



108 (x-iy 



+ 



18 (x- V 



log a; 



Z°-penguin 

C S m{x) -- 

C H {x) -- 

C£\x,y) - 

C x 2 \x,y) -. 



x (x — 6) x (3x + 2) 



--Bsm(x) 



16 (y - x) 

\fxy~ 



x y 

log x log y 



8 (y - x) 



x — 1 

X 



x — 1 



log a; 



y-i 

y 



y-i 



logy 



Gluon-penguin 

,_, x(a; 2 + 11a; — 1£ 

E S m{x) 

E H (x) 



-9a; 2 + 16a; - 4) 



12 (x - 1) 6(x-l 

a; (7a; 2 - 29a; + 16) x (3s - 2) 



log a; 



36 (x - V 



E x {x) 



-11a; 2 + 7a; - 2) 
36 (x -l) 3 



6(x-r 



x 



6(x-r 



log a; 



log a; 



17 



References 

[1] G.D. Barr et al. (NA31 Collaboration), Phys. Lett. B 317 (1993) 233. 

[2] L.K. Gibbons et al. (E731 Collaboration), Phys. Rev. Lett. 70 (1993) 1203. 

[3] A.J. Buras, M. Jamin, and M.E. Lautenbacher, Nucl. Phys. B 408 (1993) 209. 

[4] M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, Phys. Lett. B 301 (1993) 263. 

[5] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 73 (1994) 225 

[6] G. Buchalla, A.J. Buras, M.K. Harlander, M.E. Lautenbacher, and C. Salazar, Nucl. Phys. 
B 355 (1991) 305. 

[7] J.M. Gerard, W. Grimus, and A. Raychaudhuri, Phys. Lett. B 145 (1984) 400; 

J.M. Gerard, W. Grimus, A. Masiero, D.V. Nanopoulos, and A. Raychaudhuri, Nucl. Phys. 

B 253 (1985) 93; 

M. Dugan, B. Grinstein, and L. Hall, Nucl. Phys. B 255 (1985) 413; 

A. Dannenberg, L. Hall, and L. Randall, Nucl. Phys. B 271 (1986) 574. 

[8] J. Bijnens and M.B. Wise, Phys. Lett. B 137 (1984) 245; 
A.J. Buras, and J.M. Gerard, Phys. Lett. B 192 (1987) 156; 
S.R. Sharpe, Phys. Lett. B 194 (1987) 551; 
M. Lusignoli, Nucl. Phys. B 325 (1989) 33. 

[9] A.J. Buras, M. Jamin, M.E. Lautenbacher, and RH. Weisz, Nucl. Phys. B 370 (1992) 69, 
and Addendum, Nucl. Phys. B 375 (1992) 501; 

A.J. Buras, M. Jamin, and M.E. Lautenbacher, Nucl. Phys. B 400 (1993) 37 and 
Nucl. Phys. B 400 (1993) 75; 
M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, Nucl. Phys. B 415 (1994) 403. 

[10] M.J. Duncan, Nucl. Phys. B 221 (1983) 285; 

J.F. Donoghue, H.P. Nilles, and D. Wyler, Phys. Lett. B 128 (1983) 55. 

[11] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 69 (1992) 3439. 

[12] L. Ibanez and D. Lust, Nucl. Phys. B 382 (1992) ; 

V.S. Kaplunovsky and J. Louis, Phys. Lett. B 306 (1993) 269. 

[13] S. Bertolini, F.M. Borzumati, A. Masiero, and G. Ridolfi, Nucl. Phys. B 353 (1991) 591. 

[14] Particle Data Group, Review of Particle Properties, Phys. Rev. D 45 (1992) 1. 

18 



[15] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. B 209 (1988) 119. 

[16] M. Artuso et al. (CLEO Collaboration), Phys. Rev. Lett. 62 (1989) 2233. 

[17] R. Barbieri and G.F. Giudice, Phys. Lett. B 309 (1993) 86. 

[18] R. Ammar et al. (CLEO Collaboration), Phys. Rev. Lett. 71 (1993) 674. 

[19] F. Gilman and M. Wise, Phys. Lett. B 83 (1979) 83; 

B. Guberina and R.D. Peccei, Nucl. Phys. B 163 (1980) 289; 
F.J. Gilman and J.S. Hagelin, Phys. Lett. B 126 (1983) 111; 
A.J. Buras and J.M. Gerard, Phys. Lett. B 203 (1988) 272. 

[20] J.M Flynn and L. Randall, Phys. Lett. B 224 (1989) 221, Erratum B 235 (1990) 412; 
E.A. Paschos and Y.L. Wu, Mod. Phys. Lett. A 6 (1991) 93; 
M. Lusignoli, L. Maiani, G. Martinelli, and L. Reina, Nucl. Phys. B 369 (1992) 139. 

[21] G. Buchalla, A.J. Buras, and M.K. Harlander, Nucl. Phys. B 337 (1990) 313. 

[22] S. Bertolini, M. Fabbrichesi, and E. Gabrielli, Phys. Lett. B 327 (1994) 136; 
N.G. Deshpande, X.G. He, and S. Pakvasa, Phys. Lett. B 326 (1994) 307. 



19 





WW 


HH 


HW 


X 


Vi 


0.85 


0.07 


0.21 


0.62 


V2 


0.62 


0.62 


0.62 


0.62 


V3 


0.36 


0.21 


0.21 


0.62 



Table 1: Numerical values for the parameters rji which represent the QCD corrections to 
box diagrams. 



B K 


\Vcb\ 


\v ub /v cb \ 


m s (GeV) 


J l"ri+r)' 


A$ D (GeV) 


0.65 ±0.15 


0.040 ± 0.004 


0.09 ±0.04 


0.15 ±0.03 


0.25 ±0.05 


0.3 ±0.1 



Table 2: Numerical values for the relevant input parameters of the SM used in e and e'/e 



msusy (GeV) 


M H+ (GeV) 


sin 5 


n 


V, (Qe) (GeV 3 ) 


( e / e Jmin 


45 


80 


0.567 


1.08 


2.71xl0~ 2 


-6.28xl0~ 5 


45 


oo 


0.567 


0.784 


2.64xl0~ 2 


1.65xl0~ 4 


60 


oo 


0.553 


0.721 


2.69xl0~ 2 


2.10xl0~ 4 


100 


oo 


0.822 


0.745 


2.81xl0~ 2 


3.00xl0- 4 


oo 


oo 


0.978 


0.682 


2.85xl0~ 2 


4.50xl0~ 4 



Table 3: The minimum values of e'/e and the corresponding values of sin S, Q and y§ {Qe) 
for different values of m SUS y and M#+ for m t = 174GeV. The parameters are subject to the 
constraints from B° — B° mixing and b — ► 57. 



™susy (GeV) 


M H+ (GeV) 


sin 5 


n 


Ve (Qe)o (GeV 3 ) 


( e / e /max 


45 


80 


0.798 


0.446 


2.78xl0~ 2 


6.26xl0~ 4 


45 


00 


0.737 


0.355 


2.75xl0~ 2 


6.64xl0~ 4 


60 


00 


0.781 


0.435 


2.77xl0~ 2 


6.24xl0~ 4 


100 


00 


0.878 


0.555 


2.80xl0~ 2 


5.57xl0~ 4 


00 


00 


0.978 


0.682 


2.85xl0~ 2 


4.50xl0~ 4 



Table 4: The maximum values of e'/e and the corresponding values of sin S, Q and y§ {Qe) 
for different values of m SUS y and Mjj+ for m t = 174GeV. The parameters are subject to the 
constraints from B° — B° mixing and b — ► 57. 
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rasusy (GeV) 


M H+ (GeV) 


(e'/e) min x 10 4 


VMmax >< 10 4 


45 


80 


1.40± 1.30 


7.24± 2.80 


45 


oo 


3.32± 1.80 


7.30± 2.83 


60 


oo 


3.82± 2.09 


7.14± 2.80 


100 


oo 


5.00 ± 2.43 


6.83± 2.73 


oo 


oo 


6.49± 2.71 


6.49± 2.71 



Table 5: The minimum and maximum average values of e'/e and their standard deviations 
for different values of m SUS y and M#+ for m t = 150GeV. The parameters are subject to the 
constraints from B° — B° mixing and b — ► 57. 



™susy (GeV) 


M H+ (GeV) 


(e'/e) min x 10 4 


(e'AWx x 10 4 


45 


80 


-0.432± 1.04 


5.44± 2.60 


45 


00 


1.51± 1.14 


5.68± 2.66 


60 


00 


1.92± 1.33 


5.41± 2.57 


100 


00 


2.75± 1.81 


5.03± 2.46 


00 


00 


4.27± 2.29 


4.27± 2.29 



Table 6: The minimum and maximum average values of e'/e and their standard deviations 
for different values of msusy and Mh+ for m t = 174GeV. The parameters are subject to the 
constraints from B° — B° mixing and b — > S7. 



"^susy (GeV) 


M H+ (GeV) 


(e'/e) min x 10 4 


(eVe)max x 10 4 


45 


80 


-1.44± 1.28 


4.32± 2.23 


45 


00 


0.61± 0.83 


4.65± 2.43 


60 


00 


0.99 ± 0.94 


4.39± 2.32 


100 


00 


1.48 ± 1.38 


3.80 ± 2.10 


00 


00 


2.84± 1.93 


2.84 ± 1.93 



Table 7: The minimum and maximum average values of e'/e and their standard deviations 
for different values of msusy and M H + for m t = 190 GeV. The parameters are subject to the 
constraints from B° — B° mixing and b — ► 57. 
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msusy (GeV) 


M H+ (GeV) 


(e'/e) min x 10 4 


(^U >< 1Q4 


45 


80 


-0.440± 1.05 


5.39± 2.67 


45 


oo 


1.51± 1.18 


5.62± 2.76 


60 


oo 


1.92± 1.37 


5.37± 2.67 


100 


oo 


2.72± 1.81 


4.90± 2.51 


oo 


oo 


4.17± 2.33 


4.17± 2.33 



Table 8: The minimum and maximum average values of e'/e and their standard deviations 
for different values of m S usy and M H + for m t = 174GeV. The parameters are not subject to 
the constraints from B° — B° mixing and b — ■> 57. 
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Figure 1: The box diagrams AS = 2 (a-b) and AB = 2 (c-d) for e and B — B mixing 
respectively, in the minimal supersymmetric SM. 
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W,H + D (m W,s,b 



u,c,t , 



„ u,c,t 



W,H + 
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(m) 



l.Oi 



e) D 

"k. c k , t k D W 
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Figure 2: The penguin (a-d) and the box (e-i) diagrams for e'/e in the minimal supersym- 
metric SM. 
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Q 



Figure 3: The bounds of cos 5 as a function of m SU sy for the charged Higgs mass M H = 
80GeV and M H = oo respectively in figure a) and b). The parameters are not subject to the 
constraints from B° — B° mixing and b — ► 57. 
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Figure 4: The bounds of cos 5 as a function of m SUS y for the charged Higgs mass M H = 
80GeV and M H = oo respectively in figure a) and b). The parameters are subject to the 
constraints from B° — B° mixing and b — ► 57 and the solution for cos 5 < is excluded. 
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Figure 5: The bounds of e'/e with cos<5 > as a function of msusy for the charged Higgs 
mass M H = 80GeV and M H = oo respectively in figure a) and b). The parameters are 
subject to the constraints from B° — B° mixing and b — > 57 and the solution for cos 5 < is 
excluded. 
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Figure 6: The bounds of e'/e as a function of msusy for the charged Higgs mass M H = 80GeV 
and M H = oo respectively in figure a) and b). The solid and dashed curves correspond to 
the solutions with cos S > and cos S < respectively. The parameters are not subject to 
the constraints from B — B° mixing and b — > S7. 
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